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Abstract 



We introduce an adiabatic long-range quantum communication proposal based on a quantum dot array. By adiabat- 
ically varying the external gate voltage applied on the system, the quantum information encoded in the electron can 
be transported from one end dot to another. We numerically solve the schrodinger equation for a system with a given 
number of quantum dots. It is shown that this scheme is a simple and efficient protocol to coherently manipulate the 
population transfer under suitable gate pulses. The dependence of the energy gap and the transfer time on system pa- 
rameters is analyzed and shown numerically. We also investigate the adiabatic passage in a more realistic system in the 
presence of inevitable fabrication imperfections. This method provides guidance for future realizations of adiabatic 
quantum state transfer in experiments. 
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1. Introduction 

In quantum information science, quantum state trans- 
fer (QST), as the name suggests, refers to the transfer of 
an arbitrary quantum state a\Q) +/?|1) from one qubit 
to another. There are two major mechanisms for QST 
in quantum mechanics. The first approach is usually 
characterized by preparing the quantum channel with 
an always-on interaction where QST is equivalent to 
the time evolution of the quantum state in the time- 
independent Hamiltonian III, |2l [j] . However, these ap- 
proaches require precise control of distance and tim- 
ing. Any deviation may lead to significant errors. The 
other approach has paid much attention to adiabatic 
passage for coherent QST in time-evolving quantum 
systems, which is a powerful tool for manipulating a 
quantum system from an initial state to a target state. 
This method of population transfer has the important 
property of being robust against small variations of the 
Hamiltonian and the transport time, which is crucial ex- 
perimentally since the system parameters are often hard 
to control. Recently, the adiabatic method has been 
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applied to a variety of physical systems to realize co- 
herent QST. Among these, the typical scheme for co- 
herently spatial population transfer has been indepen- 
dently proposed for neutral atoms in optical traps |01 
and for electrons in quantum dot (QD) systems 15J via a 
dark state of the system, which is termed coherent tun- 
neling via adiabatic passage (CTAP) following |5J]. In 
such a scheme, the tunneling interaction between adja- 
cent quantum units is dynamically tuned by changing 
either the distance or the height of the neighboring po- 
tential wells following a counterintuitive scheme, which 
is a solid-state analog of the well-known stimulated Ra- 
man adiabatic passage (STIRAP) protocol 16[] of quan- 
tum optics. Since then, the CTAP technique has been 
proposed in a variety of physical systems for transport- 



, electrons 
HQ. It has also 



ansport 



ing single atoms J7|,|8|], spin states ||9 
and Bose-Einstein condensates Jl 
been proposed as a crucial element in the scale up to 
large quantum processors IU5[|l6ll . 



Recently, Ref. [11711 presented a scheme to adiabati- 
cally transfer an electron from the left end to the right 
end of a three dot chain using the ground state of the 
system. This technique is a copy of the frequency chirp- 
ing method 11181 Il9ll . which is used in quantum op- 
tics to transfer the population of a three-level atom of 
the Lambda configuration. The scheme Il7ll is pre- 



Preprint submitted to Sci China Ser G-Phys Mech Astron 



June 29, 2012 



O^-o o 



(')# 



Figure 1 : Schematic illustrations of adiabatic QST in a multi-dot ar- 
ray. The system is controlled by gate voltages, /i a (t) (a = A, B). By 
adiabatically varying the gate voltages, one can achieve long-range 
QST from the left end to the right QD of the proposal. 



sented as an alternative to a well known transfer scheme 
(CTAP) H. However, different from the CTAP process, 
the protocol in Ref. Il7ll considers a three QD array with 
an always-on interaction that can be manipulated by the 
external gate voltage applied on the two external dots 
(sender and receiver). Through maintaining the system 
in the ground state, it shows that it is a high-fidelity 
process for a proper choice of system parameters and 
also robust against experimental parameter variations. 
In this paper we will consider the passage through the 
Af-site coupled QDs array (tight-binding model), which 
is schematically illustrated in Fig. Q] Gates applied on 
the two end dots control the on-site energy of each dot. 
In particular, the nearest-neighbor hopping amplitudes 
are set to be uniform. We first investigate the effect 
of system parameters on the minimum energy gap be- 
tween the ground state and the first excited state. Tak- 
ing a 5-dot structure as an example, we show that the 
electron can be robustly transported from one end of the 
chain to the other by slowly varying the gate voltages. 
This structure is easy to extend to an arbitrary number 
of sites. 

The paper is organized as follows. In Sec. II the 
model is setup and we describe the adiabatic transfer of 
an electron between QDs. In Sec. Ill we show numer- 
ical results that substantiate the analytical results. The 
last section is the summary and discussion of the paper. 



2. Model setup 

We introduce a simple tight-binding chain with uni- 
form nearest-neighbor hopping integral -J as a quan- 
tum data bus for long-range quantum transport, see 
Fig. Q] The sender (Alice) and the receiver (Bob) can 
only control the external gate voltages /i ff (?) {a = A, B), 
which are applied on the two end QDs. In this proposal, 
the quantum information cos 6 |f) + sin 9 1|) encoded in 
the polarization of the electron can be transported from 



Alice to Bob. In this scheme the spin state of the elec- 
tron is a conserved quantity for the Hamiltonian of the 
medium, so the spin state cannot be influenced during 
the propagation. For simplicity, we will disregard the 
electron's spin degrees of freedom in the following dis- 
cussion and just illustrate the principles of QST. 

Accounting only for the occupation of the lowest 
single-particle state of each dot, the system is described 
by an A^-site tight-binding chain. The Hamiltonian can 
be written as 



'Hu = { a ) a i +i + hx -) ' 

<H C = Ha (0 a\a\ + [x B (0 a\a N , 



(1) 



where <r(a ; ) denotes the spinless fermion creation (an- 
nihilation) operator at the j'-th quantum site. -J (J > 
0)is the coupling constant, which accounts for the hop- 
ping of the electron between dots j and j + 1. The 
on-site energies /44(f) and /ig(f) are modulated in Gaus- 
sian pulses to realize the adiabatic transfer, according to 
[shown in Fig. 2(a)] 

fi A (t) = -HA,max exp [-a 2 ? 2 ^] , (2a) 
(1b(0 = -HB.max exp [-a 2 (t - t) 2 /2] , (2b) 

where r and a are the total adiabatic evolution time and 
standard deviation of the control pulse. For simplicity 
we set the two peak values to be equal, i.e., fiA,max = 
HB,max = Ha (m > 0), in the discussion that follows and 
choose //() » /. 

In this proposal, we will concentrate on the single- 
particle problem and use the ground state |^ g (/)) of the 
Hamiltonian H(t) to induce a population transfer from 
state 1 1 ) to |A0. Starting from t = 0, we have //4(C)) = 
-//o and /<b(0) ~ 0. The Hamiltonian at t — reads 

N-l 

<H(f = 0) = -//oajfli - J ^ (a\aj+\ + h.c.) . (3) 
j=i 

The ground state of Eq. © is a bound state, which can 
be obtained via the Bethe ansatz method. A straightfor- 
ward calculation shows that 



(4) 
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where £ = 7///o- By choosing a sufficiently large 
value of fio, the ground state 1^(0)^ can be reduced to 
1^,(0)^ * |1) = a[ |0). To illustrate with an example, 
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the probability of |1) in \i// g (t = 0)^ can achieve 99.75% 
when the peak voltage is set to be fio/J = 20. 

With the same reasoning, in the time limit t = t, the 
parameter //a(0 goes to zero and /Ub(0 goes to -fiQ, Due 
to the reflection symmetry (relabelling sites from right 
to left) of the system, we can see that 

/ N 

l*,(f=T)) = yji-e^ N - j \j). (5) 

One can see that the ground state of Eq. ([1]) evolves 
to be \4>g{t = t)\ ~ \N). Preparing the system in state 
14* (f = 0)) = |l)and adiabatially changing //a(0 and 
IJ-nit), one can see that the system will end up in \N) 

\W(t = 0)) = \l)^\W(t = r)) = \N}- (6) 
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Figure 2: (a) Gate voltages (in units of /io) as a function of time (in 
units of r) described in Eq. (2). >m(0 is the solid line and /ig(f) is 
the dashed line, (b) The instantaneous eigenenergy of the lowest two 
states ip\ and if/ g though the gate pulse shown in (a) for the values of 
fJo = 20, / = 1.0, and a = 5/rinan/V = 5 structure. The gap is 
minimum at / = t/2, A = £i(t/2) - e s (r/2). 

The analysis above is based on the assumption that 
the adiabaticity is satisfied. The crucial requirement for 
adiabatic evolution is 

\s g (t) - El (t)\ » 1(^(01^1(0)1 , (7) 



which greatly suppresses the quantum transition from 
the ground state \i/r g (t)) to the first-excited state \iffi(f)). 
Firstly, one must make sure that no level crossings oc- 
cur, i.e., s g (t) - ei(0 < 0. To evaluate instantaneous 
eigenvalues of the Hamiltonian is generally only pos- 
sible numerically. In Fig. |2{b) we present the results 
showing the eigenenergy gap between the instantaneous 
first-excited state and ground state undergoing evolu- 
tion due to modulation of the gate voltage according the 
pulses given in Eq. (2) for fig = 20, 7=1.0 and a = 5/t. 
The eigenvalues shown in this figure exhibit pronounced 
avoided crossing and approach nonzero minimum A = 
£i(t/2) - s s (t/2) at t — t/2. This minimum energy gap 
plays a significant role in the transfer, because the to- 
tal evolution time t should be large compared to 1/A. 
In this scheme, the energy gap A depends both on the 
number of QDs and gate voltages. To study the rela- 
tionship between the total evolution time r and system 
parameters is one of the important contributions of this 
paper. 

Fig.[3]shows the effect four factors have on the energy 
gap A. In Fig. a), for a given a — 5/t and J = 1.0, we 
plot the energy gap A as a function of N~ 2 for jUq =16, 
20, and 24. The numerical results indicate that the gap 
A ~ J 2 / (fioN 2> ), which implies that the adiabatic trans- 
fer time of QST grows quadratically with the spatial 
separation of the two end states because the minimum 
gap plays an opposite role for the adiabatic QST. The 
other thing is that A is also determined by the dimen- 
sionless parameter a. As an example, Fig. [3jb) shows 
the numerically computed behavior of A as a function of 
a with N - 5 and /io// = 20. From Fig. [3b) we see that 
the gap becomes larger as a increases and then tends to 
be a constant 0.732, which is the gap of tight-binding 
chain (Af = 5) without any on-site energy. The reason is 
that the bigger <rr is, the smaller the overlap amplitude 
of the two pulses, i.e., exp(-a 2 T 2 /8) — > 0. The energy 
gap A of t H(j/2) then approaches the maximum value 
2J[costt/(N + 1) - cos 2tt/(N+1)]. 

3. Numerical Examples 

In this section let us firstly review the transfer process 
of this protocol. At t — we initialize the device so that 
the electron occupies site-1, i.e., the total initial state is 
\V (0)) = 1 1 ), and slowly apply gate pulses, which re- 
sults in robust transport of the electron from one end of 
the chain to the other. The consequent time evolution of 
the state is given by the Schrodinger equation (assuming 
h = 1) 

4pP(0> = 7f(0rF(0>. (8) 
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Figure 3: Effect the system parameters have on the energy difference, 
(a) The gap A obtained using a numerical method for the systems 
N = 5, 6, 7, 8, 9, and 10, with /i = 16, 20, 24 and J = 1.0 are plotted. 
It indicates that A ~ J 1 1 (ji N 2 ). (b) The gap A as a function of a. It 
shows the increase of the gap with increasing a. 



The time evolution creates a coherent superposition: 



\y(t)) = Y j C j (t)\j): 



(9) 
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where Cj(t) denotes the time-dependent probability am- 
plitude for the electron to be in the j-th QD that obeys 
the normalization condition "Z/J =1 \cj(t)\ = 1. At time t 
the fidelity of the initial state transferring to the dot-N is 
defined as 



F(r) = |<AWr)>| 2 = Mr)| 2 



(10) 



A feasible proposal should be able to perform effi- 
cient high-fidelity QST in the shortest possible time. In 
order to provide the most economical choice of param- 
eters for reaching high transfer efficiency, we used stan- 
dard numerical methods to integrate the Schrodinger 
equation for probability amplitudes. Because the 



scheme relies on maintaining adiabatic conditions, we 
examine the effect of system parameters on the target 
state population. 

We show in Fig. [4] the probabilities as a function of 
time for different values of a where we take a 5-dot 
structure for example. The time behavior of /14(f) and 
jUs(f) follow Gaussian functions with {Iq/J = 20 [see 
Fig. 2(a)] and have been performed in a finite time. 

For the a = 4/r case with r = 500, which departs 
from the adiabatic limit, we find the result in Fig.[4ja). 
Here the population transferred to the target state is only 
about 11%. The reason is that the population is excited 
to the upper energy states through nonadiabaticity. It is 
necessary to point out that if one enlarged t extremely 
in this case, adiabaticity would be fulfilled, which would 
result in high transfer fidelity. 

On the other hand, as shown in Fig. [3] enlarging a can 
increase the level spacing between the first excited state 
and the ground state and hence cause the adiabaticity of 
the system to become better. In Fig. |4|b) we show the 
population evolution taking a — 5 /t as an example and 
the populations of the states |1) and \N) are exchanged 
with a fidelity of 99.5%. 

The choice of pulse modulation is therefore impor- 
tant with the maximum transfer speed ultimately con- 
trolled by the adiabatic criteria for the transfer. To see 
the emergence of the adiabatic limit, we plot in Fig. 4(c) 
the transfer fidelity as a function of the adiabaticity pa- 
rameter a. One can see that as a increases there is an 
exponential appearance of the adiabaticity in the ideal 
limit. That means the smaller the overlapping of two 
pulses is, the more ideal adiabatic transfer takes place. 
This result is extendible to an arbitrary number of QDs. 
In the discussion that follows, we choose a - 5/t. 

In section II, we showed that the energy gap A also 
depends on the system parameters, such as the number 
of QDs N and coupling strength J. In order to quantita- 
tively determine the time needed to achieve high fidelity 
QST, we solve the schrodinger equation for constants 
//o = 20, J — 1.0. In Fig. [5] we present results showing t 
as a function of the QDs number N. The quadratic curve 
fitting shows that the minimum possible transfer times 
are proportional to A^ 2 , giving a high-fidelity transfer of 
\c?j(t)\ 2 > 0.995. On the other hand, the energy gap 
A decreases when the peak value rises [see Fig.[3ja)]. 
Consequently, the time scale t is proportional to and of 
the order of fio/J 2 for a given N. To sum up the above 
discussion, in practice the minimum possible transfer 
timescale of this adiabatic passage will be of the same 
order as N 2 /io/J 2 . 

For a long enough evolution time r, the maximum 
fidelity of this scheme depends on the contrast ratio be- 
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Figure 4: Probability of rinding a single electron in basis states in the 
ground state as a function of time for the values of N = 5, J = 1.0, 
Mo = 20, r = 500 and (a) a = 4/r, (b) a = 5/t. (c) Transfer fidelity 
|c;v(f)| 2 as a function of parameter a. As o- is increased it is more able 
to obtain high fidelity transfer. A smaller a introduces a population of 
excited states and the transfer is no longer complete. 
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Figure 5: Total evolution time t (in units of fio/J') as N is increased 
under the condition that \cn(t)\ 2 > 0.995. As Wis increased the energy 
gap is decreased, resulting in longer evolution time across the array. 
The solid line is the quadratic curve fitting, which indicates that the 
evolution time grows quadratically with the number of QDs. 



tween peak values n Kimax (k = A, B) and coupling con- 
stants J. The reason is that small peak values improve 
adiabaticity, but lead to a low fidelity because the initial 
and final energy eigenstates are not the desired states |1) 
and \N), respectively. 

To determine the parameter range of \i Ktmax (k = A, B) 
needed to achieve high fidelity transfer, we numerically 
integrate the density matrix equations of motion, with 
varying peak values /J. A ,max/J and Hn.max/J from 10 to 
25. Fig. [6] shows the transfer fidelity \cn(t)\ 2 plots as 

a function Of /JA,max/ J and ^B,max 

I J for t — 1000 and 
a — 5 /Yin a 5-dot system. The fidelity approaches unity 
as fiA,max and HB,max increase. The figure is nearly sym- 
metric with respect to the line ^A,max = HB,max- Fig- IS a) 
is taken from Fig. |Sb) by slicing through the diagonal 
gray line. We can see that to realize near-perfect fidelity 
transfer (\c n (t)\ z > 0.995) one has to use peak values 
satisfying pL KMax > 207 (a: = A, B). 

The other advantage of this scheme is defect toler- 
ance of the system parameters. We now assume that the 
tunnel coupling has a random but constant offset Ssj, 
i.e. Jj = /(l - Ssj), where sj is drawn from the stan- 
dard uniform distribution on the open interval (0, 1) and 
all Sj are completely uncorrected for all sites along the 
chain. We show some examples in Fig. [7] for QST in 
the chain of N = 5 with maximum coupling offset bias 
6 = 0.1,0.2, and 0.3. It shows that weak fluctuations 
(up to 6 = 0.2) in the coupling strengths do not deterio- 
rate the performance of our scheme. For 6 — 0.3 we can 
see that arbitrarily perfect transfer remains possible ex- 
cept for some rare realizations of 0.3s. To realize high- 
fidelity QST transfer, the price for unprecise couplings 
is thus a longer transmission time. 
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4. Summary 




Figure 6: (Color online) Plot of the transfer fidelity \cn(t)\ 2 as a func- 
tion of peak values (in units of 7). (a) HA/nax = HB,max = j"o varies 
from 107 to 257; (b) iiajmu and UBjnax vary from 107 to 257 respec- 
tively. The contour lines, labeled with the corresponding values of 
|cjv(t)| 2 , display the gradual increase of transfer fidelity as jiA,maxlJ 
and HBjnaxlJ grow. The fidelity is close to one (kiv(r)| 2 > 0.995) 
when two peak values are achieved for fi Ktmax > 207 (k = A, B). 
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Figure 7: (Color online) The transfer fidelity |cw(r)| 2 for a tight- 
binding chain of length N = 5 with ytAmax — HB.max = 20 J, t = 500 
and a = 5/r. The coupling strengths are chosen randomly from the 
interval [(1 - <5)7, 7] for <5 = 0.1 (Square), S = 0.2 (Circle) and <5 = 0.3 
(Triangle). The number of random samples is 20. 



We have introduced a robust and coherent method of 
long-range coherent QST through a tight-binding chain 
by adiabatic passage. This scheme is realized by modu- 
lation of gate voltages applied on the two end QDs. Un- 
der suitable gate pulses, the electron can be transported 
from one end of the chain to the other, carrying along 
with it the quantum information encoded in its spin. 
Different from the CTAPn Scheme l2(ill . our method is 
to induce population transfer through the tight-binding 
chain by maintaining the system in its ground state and 
this is more operable in experiments. We have stud- 
ied the adiabatic QST through the system by theoretical 
analysis and numerical simulations of the ground state 
evolution of the tight-binding model. The result shows 
that it is an efficient high-fidelity process (> 99.5%) for 
a proper choice of standard deviation a > 5/r and peak 
values yUo > 20 of gate voltages. For an increasing num- 
ber of dots, we found that the evolution time scale is 
r ~ N 2 fio/J 2 - We also consider the QST along the quan- 
tum chain if their coupling is changed by some random 
amount. We further find that weak fluctuations in the 
coupling strength still allow high fidelity QST. 
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